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Examining an approximation inspired by equalization theory, we 
consider the minimum angle SIn between the subspaces ofHilbert space 
generated by the sequences \e' k,j: \ N -N and [e ika }\k\ >jv« Here a> £ 
[— 7r,7r] and the inner product for the Hilbert space involve a positive, 
bounded weight function r(cc). The finite Toeplitz matrices R and T 
generated by r(co) and l/r(u)), respectively, play a crucial role, and, in 
fact, sin 2 ttN is the reciprocal of the largest eigenvalue ofRT. In general, 
sin 2 QN is shown to be bounded away from unity as N becomes large. 
The geometry of the problem enables us to give some results concerning 
the product matrix RT, which, out of the present context, may seem 
surprising. 

I. INTRODUCTION AND SUMMARY 

Let H be the Hilbert space of square-integrable functions on [— ir,ir] 
with an inner product given by 

(f,g)r = ^- f*/*(«)*(co)r(co)dw, (1) 

where the weight function r(o>) is bounded and strictly positive; i.e., 

0<r±r(a>)^R. (2) 

We call a Fourier subspace of H any subspace generated by a finite 
or infinite collection of functions of the form e inu , n an integer. In par- 
ticular, we shall be interested in the Fourier subspaces [relative to the 
metric r(a>)]: 

(3) 
G N = I E g n e ifm 

\\n\>N 



tThe subscript r, as in (1), will be used when we wish to emphasize that the weight 
function r(w) is being used. No subscript will refer to an arbitrary inner product, while 
the case r(w) = 2ir will be called the "usual" metric. The usual inner product will be written 
without a subscript as well. 
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for N ^ O.t If r(w) is constant, then the subspaces F N and G N are or- 
thogonal. We will be concerned with the minimum angle between F N 
and Gn for a general weight function satisfying (2), and with the limiting 
behavior of this angle as N -* «. (The concept of the angle between 
subspaces is not new to the engineering literature. See for example Ref. 

1-) 

The main results of our investigation are stated in terms of two finite 

Toeplitz matrices, R and I\ which are generated by the weight functions 

r(w) and g(a) = l/r(«), respectively [see eqs. (30) and (33) for precise 

definitions]. We also need the Fourier coefficients r„ and g„ of r(a>) and 

g(w). Then, we show 



(£) sin 2 Q n = 
( 



largest eigenvalue of R r 

[1 /•* 1 (** do) ~| — i 

1+— ( r(co)dcoX— I — " <1. 

(Hi) All eigenvalues of RT are ^ 1. 

0») f |»|r„^;<0. 

This is not a trivial inequality in the sense that £ -o,r n ^„ = 1 is. 

The case r(w) = 1 + a cos a> is solved exactly in Section V, showing that 
the obvious bound sin 2 fiA/ s r min /r max is often loose. A better bound, still 
involving only this ratio, is given in (68). 

In somewhat general terms, this problem arose in the mean-square 
equalization theory of data transmission, where the question is one of 
bounding the effect of replacing tap weight values by certain Fourier 
coefficients. To be specific, let us ignore the effects of noise and note that 
the job of the equalizer is to invert the Nyquist equivalent channel. That 
is, if we had an infinite number of taps at our disposal, we would take the 
transfer function of the equalizer C m (a>) to be 

C(«)--=7-r- £«nC inM . 

A(co) -« 

The equalizer transfer function C N ( ) when only the usual (22V + 1) taps 
are available can always be written as 

Cn(«)--~ t-EM'-- E m* 1- . 

X(w) -N \n\>N 

In the above expression b„, \n\ < N are "corrections" to the Fourier 
coefficients t n , \n\ < N. The mean-square error resulting from the 

t Any function of the form (3) is in H if and only if the associated coefficient sequence is 
square-summable. 
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2 |X(co)| 2 da;. 



equalizer CW(w) can then be shown to be given by 

7-£|X(«)c N («)-i|»a« 

1 /•*■ I N 

The minimum mean-square error E 2 min is the minimum of the above 
expression over the 5/,. Now we can imagine taking the (fixed) vector 

£ *ke iku 

\k\>N 

and decomposing into a vector in the space Fyv and one perpendicular 
to it. The part in F N can be "subtracted off by the choice of 5's, leaving 
the remainder. The fraction subtracted off can never be greater than 
cos 2 fi/v, where tt N is the angle between the two subspaces F N and G/v 
when X eq (o;) 2 is used as the weight function r(o>) for inner products. 
Thus, 

sin 2 fi N X || e N || 2 < E 2 min ^\\e N \\~, (4) 

where 

2ir «/-7r | \k\>N 

The point of replacing the exact tap values of the finite equalizer by 
Fourier coefficients is not to replace one calculation by another. Rather, 
it is to supplement calculation by insight, since much is known about 
properties of Fourier coefficients. This will be done for a specific equa- 
lization problem involving timing recovery for finite equalizers in a later 
work. 

II. WARMUP EXERCISES 

The angle 6 between two fixed vectors, / and g, is defined by 

9 = C0S "'wS' * e[M - ,5) 

If f and g are restricted to be in subspaces F and G, respectively, the 
infimum of (5) (call it ft) over all / and g (so restricted) is called the angle 
between the two subspaces. We easily see that 

11/ -g\\ 2 ^ 11/11 2 + Wl 2 - 2||/|| X ||*|| cos Q, (6) 

and thus by minimizing the right member of (6) with respect to the norm 
of /, we have 

inf ll/-£|| 2 i=sin 2 n||s|| 2 . (7) 
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In fact, we can also calculate sin 2 fi via the formula 

sin 2 fi = infinf — 5— ;rr — . (8) 

8 ( llS II 2 

When (2) holds, the infimum angle between our subspaces F N and Gn 
[given by (3)] is actually attained and its value is strictly positive. In fact, 
it follows from an application of a theorem by Paley and Wiener 2 that 
the two sequences (usual metric) 



fonl 88 



vA^H 



(9) 



2x 

v27rr(oj) J 
form a complete biorthogonal pair, i.e., 

(4>n,fm) = Snm (10) 

and 

k - £ (*„,/»)*„ = £ (<An»0n (ID 

for any /i in L 2 . Thus, either sequence in (9) forms a basis for L 2 , or, 
equivalents, {e' nw } r forms a basis for H [with weight function r(»)]. Now 
iffEF N ,gSG N , 

inf|l/-g|lr 

is attained when g is the orthogonal projection of/ on Gyv, and is a con- 
tinuous function of the finite dimensional /. Therefore, 

illf mf ||r||2 ' 

( L g 11/11; J 

is attained, since we may restrict ||/|| r = 1 and thus are minimizing over 
a compact set. The basis property of \e inu \ r in H assures that the mini- 
mum is not zero. 

There are several ways to get at the minimum angle Un between Fn 
and G N . We shall begin by using (8) and the calculus of variations. 

However, before we begin to work on this, let us review an old problem 
of linear prediction (really, linear interpolation) theory. We are required 
to find the minimum value of 

E 2 = — f*|l- I ame' m 1 r(to)d<j (12) 

over all h sequences \a m ), under assumption (2) for r(w). We let 

a(a>)= E a m e im » 
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be any element of Li{— ir,ir) which has zero for its zeroth Fourier coeffi- 
cient. We then have 



E 2 = ^- f*|l-a(tt)| a r(«)d«. 



(13) 



The calculus of variations yields 

f8a*[l - a(a)]r(<a)da = (14) 

or 

fe in »[l - a(a>)]r(co)da; =0, n ^ 0. (15) 

Thus, 

[1 - a(w)]r(«) = const = k. (16) 

From (13) and (16), on the one hand, we have 

„ 2 \ r* k 2 k 2 /•» du 

and, on the other, 

Emin = — f* d«(l - a*)[(l - a)r] = £- f* (1 - a*)da> « ft, (18) 

■&7T •/— ir Z7T «/ — k 

since a* has no m = term. Equating the results of (17) and (18) enables 
us to solve for k [note that k = must be excluded under (2)], yield- 
ing 

£min = — =k. (19) 



1 /•*■ do) 

2rr J-* r(oj 



•(o>) 

If r(co) has a zero somewhere in such a manner that fl/r = «>, then (19) 
says Emi n = 0. This turns out to be the correct conclusion for the infimum, 
but this infimum is never attained. As is well known, the calculus of 
variations can only be applied if the infimum is attained. In fact, in the 
present problem, if we set k = in (16), we would conclude that there 
is an li sequence \a m \ such that 

1 - E a m e imu = a.e., (20) 

which obviously cannot be. 

Another way to do this problem is to use the biorthogonal sequences 
(9). We note that all vectors of the type 

£ a m <t> m (21) 
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form the subspace orthogonal to the vector ^ . Hence, 2?m in must simply 
be the squared norm of the projection of 0o onto \po. This is (in the usual 
norm) 

IU 112 2^a f\- IU 112 I (0O>^O) 1 2 
||0OI|COS Z (0O»V'O) - II 00 II* II . II o II / ||2 

1 00 II II WO II 



/— f* Vrfa)x-==d«Y 
\2ir J-* Vrja) I 



1 /»"■ da 1 p* rf^ 

2tt J-* r(«) 2tt J-tt r(w) 



(22) 



III. FINDING THE MINIMUM ANGLE 

We proceed with the calculus-of- variations approach to finding sin 2 Q;v 
via (8). Let 

f(a)=Zfke ik »eF N 

-N 

g(a) e G N . (23) 

Then, if we vary g* in \\f - g \\ 2 „ we obtain 

J 8g*(w)[f(u) - g(a)]r(a)da = (24) 

for all allowed variations. Thus, (24) means \f - g]r E F N , or, in other 
words, 

\f((a) - g(a)]r(a) = L b k e iku = b(a) ' (25) 

-N 

for some numbers bk- As in Section II, (25) permits us to write two ex- 
pressions for min ||/ — g || 2 r . They are 



and 



. ... ||2 1 f* |M*>)l 2 , 
mm\\f-g\\ 2 r = — I - ' da 

g 27T J-ir r(oj) 

min|/-g|;-^- r X dco/*(o))6(a;) = f*-b. 



(26) 



(27) 



The vector notation in (27) refers to a row of (22V + 1) numbers. Letting 
b = k(f + bj_), f* • b j. = 0, we may equate (26) and (27), solve for k, and 
obtain 

/ N \2 
I If _,||2 • 

mm ||/ -g\\r= N N 

Zfke ikm + E Y*e*'*" 



min— C —^ — *L da (28) 

y 2% %)-* r(a) 



y 

f-7 = 
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Thus, from (8) and (28) 



sin 2 ft/v = 



mm 

u 



N 
E|Ui|2«l 

-N 



1 /»jt I N 2 I stir 

— I E u^e'*" r(«)d« X min — I dw 

27T «/— 7r I -/V w 27T «/-ir 



iV 

E (u A + w*)e'*« 

-A/ 



r(«) 



(29) 



If we introduce the Toeplitz matrix 

r nm = ^- C eiin ~ m) "^\ do3 > \*l \™\*N, (30) 

the second term in the denominator of (29) has the form 

u+Tu + w + Tw + 2u + Tw. (31) 

Expression (31) may be minimized over the appropriate w using a 
Lagrange multiplier, yielding 



u+T-tu 



Thus, 



. 2o . u+r- 1 !* 

sin z i2yv = mm 



(32) 



u+Ru ' 
where R is the Toeplitz matrix corresponding to r(co), i.e., 

Rnm = -^- C" e^-^r^da, \n\, \m\ ^ N. (33) 

2.TT %) — Jr 

Finally, the minimization of (32) yields 

Theorem I. Let the matrices R and T be defined as in (30) and (33). 
Then the minimum angle Q^ between the Fourier subspaces F^ and 
Gn, defined in (3), is 



sin 2 Qw = 



(34) 



largest eigenvalue of R F 

This theorem implies that Qjv is invariant under the replacement r(co) 
- l/r(a>). 
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Since similarity transformations preserve eigenvalues, we note that 
the eigenvalues of R T are the same as those of Vf R Vf, which is, by (2), 
(30), and (33), a strictly positive definite Hermitian matrix. 

Before exploring consequences of (34), we shall rederive it from a 
geometric point of view. Let [0/}f and l^,)r be complete biorthogonal 
sequences of vectors for a Hilbert 'space. Let 

V-lfclf, W=\<t>i\N + i (35) 

be subspaces generated by the indicated vectors. Note that the orthog- 
onal complement of W, W L , is given by W ± = T. Also note that our 
problem is equivalent to that of finding the minimum angle between V 
and W. If v £ V, and a is the angle between v and W (i.e., the angle be- 
tween v and its projection on W), and /3 is the angle between v and W ± , 
we have t 

a + = f' (36) 

Thus, the minimum angle between V and VK (call it fi) is the complement 
of the maximum angle between V and T, called 0m- Thus, we have 

sin 2 fi = cos 2 0a/. (37) 

The spaces V and T both have dimension N here. 

Let P represent the orthogonal projection operator onto T, and Q the 
orthogonal projection operator onto V. It can be shown that if V £ V 
is a vector in V which attains the minimum or maximum angle between 

V and T, we must have ' 

QPu = \v, (38) 

where X ^ is the square of the cosine of the indicated angle. We shall 
see shortly that (34) is a form of (38) when we represent P and Q by 
matrices that are representations of the restrictions of P to V and Q to 

T. 

We begin by deriving these matrix representations. For general 
biorthogonal sequences |0,j, {&■}, let 

Rnk = (<t>n,<t>k) n>k = h2 ,...,N. (39) 

Tnk = Wn,M 

t If we call the projection of u on W by ui, then v - Vi is the projection of v on W ± . Since 
v, vu and u - u, all lie in a plane, (36) follows immediately from a simple diagram depicting 
these three vectors. 

' Let a be any vector in the space such that Qa ^ 0. Then if 6 is the angle between Qa £ 

V and W, we have cos 2 = \\PQa\\ 2 /\\Qa\\ 2 = (a,QPQa)/(a,Qa). Vectors « which yield 
stationary values of this ratio of quadratic forms can be obtained by differentiating 
(a,QPQa) holding (a,Qa) constant (via a Lagrange multiplier X). This procedure yields 
(38) upon setting Qa = v. 
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Any vector x can be written uniquely as a vector in V plus a vector in the 
orthogonal complement of V, i.e., 

x = E Oi0i + f Mi. (40) 

1 iV+l 

If we form the inner product of (40) with </>;,;' = 1,2, • • • , N, we can cal- 
culate 

a* = E (R~ l )ki(<t>i,x). (41) 

/=i 

Thus, given any x, its projection onto V is simply 

N 

T.ai4>i (42) 

l 

with a, given by (41). Similarly, the projection of x onto T is 

E Mi (43) 

i 
with 

&A = E (r-^w^x) (44) 

/=1 

Hence, if we start with any vector u £ V, 

N 

^ = E Ui0;, (45) 

l 

the result of projecting it onto T and then projecting this vector back 
to V is another vector v" £ V with components v" m given by 

i>« - Z CR- 1 r- 1 )mi»i. (46) 

Hence, the operator equation (38) becomes the N X N matrix equa- 
tion 

(r#) _1 u = Au- (47) 

Equation (34) is thus rederived, after an appropriate relabeling of in- 
dices. 

Since the reciprocals of the largest and smallest eigenvalues of RT have 
interpretations as squares of cosines of angles, we have 

Theorem II. Let matrices R and T be defined as in (30), (33). Then all 
eigenvalues of RT are si. 

IV. IMPLICATIONS CONCERNING SIN 2 ft N 

From (30) and (33), we see that the matrix elements of rT are given 
in terms of the Fourier coefficients r, and gj of r(a>) and g(u) = l/r(oj). 
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Thus, (33) reads R nm = r m -n. We see that, in this notation 



N 



(RT) nk = £ gk- m r m -n,\nl\k\*N. (48) 

m = -N 



This can also be written 



• 2AT+1 , ,, 

sin — - — (co— a) ) 



(RD„*=- 1 - f X e «—*»',lM ? dalda/ . (49) 

(2 ttYJ-t r{a> ) , 1 

sin - (a?— a> ) 

Equation (49) follows easily from (36), (33), and the identity 

. 2N+1 . 
sin — - — y 

N ■ , 2 



m = -N . \p 

sin r 

2 



If the largest eigenvalue of RT approaches 1 as N becomes large, then 
sin 2 fi,v — 1, and the subspaces F N and G N described in (3) eventually 
become orthogonal. Equivalent^, we have seen that the question be- 
comes the following: Does the largest angle between the subspaces (with 
the usual metric) 






....... N 

F N = 

■N 



and 

1 N 



Gh - 



inw 



V2irr(oj) 



-N 



approach zero? If we set N = °°, the two generated spaces are identical 1, 
(all of L 2 ), so from this point of view it comes as a surprise that the lim- 
iting angle between Fn and Gyv is bounded away from zero. 

Le,t us assume that r(a>) has only a finite number of Fourier coeffi- 
cients; that is, assume rj ■ if \j\ > k. For this case, the reader may 
verify, using either (48) or (49), that the (2N + 1) X (22V + 1) matrix fir 
has the form (once N ^ k)* 



t F N is never G N for finite N unless r(u) = const. This follows from using (11) to show 
that you cannot expand each <t>„, n ^ |AT| in terms of the \pk, \k\ < N. 

t To see this, let us evaluate (flr) ab from (48). If \a | < (JV - k), the summation in (48) 
may be extended from JV = -»to7V = +», since rj = if i > k. Using the duality between 
hi and L 2 , the resulting sum is then 1/2tt fl r r(w)exp(iau))][g(w)exp(ibw)\* dw = b ab , since 
g(w) = l/r(w). 
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A X X X 

1 

1 



O 



o 



c x x x 



x 


B 


1 




X 


D_ 



(50) 



That is, the first k rows and the last k rows are nonvanishing. The 
remaining 2{N - k) + 1 diagonal elements are exactly unity, while all 
other matrix elements vanish. The four k X k matrices in the corners are 
singled out for special attention and are labeled A, B, C, D. The capital 
X's in the first and last rows are inserted only to indicate that elements 
in the first and last k rows are not vanishing, in general. 

As an example, we write the elements of A explicitly, labeling the el- 
ements of A by a rs , r, s = 0,1, • • • , k - 1; i.e., aoo = (RY)-n - N , etc. 
Then 



£ fie 



iluin—i 



.— C* l= ~ r 

^ " 27T J-n /(u>) 



da). 



(51) 



The elements of D can be determined from A using the property 

(RTy mi i m (RT)- m ,-i. (52) 

It is important to note that (for N ^ k) the elements of A and D do not 
depend on N. However those of B and C do. For example, the upper right 
corner of B is the element (Rr)-M,N, given from (48) as 



(RF)- N ,n - E fig 2 N-l- 

1=0 



(53) 



Not only does (53) depend on N, but, by the Riemann-Lebesgue lemma, 
S2N-1 ~ *" as N increases for / bounded. As N increases, all elements of 
B and C similarly vanish. 

We now look further at the problem of calculating the eigenvalues of 
(50). If one is not an eigenvalue then the matrix RF — XI has 2(n — k) + 
1 diagonal elements (1 — X). By multiplying a row in which such an ele- 
ment occurs by the appropriate constant and adding the result to one 
of the first or last rows, all the elements indicated by "X" in (50) can be 
made to vanish. Clearly then, the eigenvalues that are not unity are given 
by the nonunity eigenvalues of the 2k X 2k matrix 



K 



J A B ] 
lc DY 



(54) 
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where B and C depend on N. All other eigenvalues of (49) are one. Since 
B and C vanish in the limit of large N, we .Jiavet 

lim X m ax(N) = largest eigenvalue of A. (55) 

TV— oo 

The simple fact that F N is never Gn for finite W implies that the largest 
eigenvalue of RT, and hence K, is strictly greater than one. Also the fact 
that all eigenvalues of fir are greater than or etjual to 1 implies tr K = 
tr A + tr D > 2k. But A and D do not depend on\iV and have the same 
trace. Hence, tr A > k, and A has an eigenvalue\trictly greater than 
unity. Thus, from (34) and (55), lim sin 2 fi;v < 1. 
The above discussion implies the following:* 

Theorem III. Ifr(o)) has only a finite number of nonvanishing Fourier 
components and is not constant, then Zirajv— » sin 2 Q.N < 1- 

Theorem IV. Let 1 ^ r(w) > on [— w,ir] have only a finite number of 
nonvanishing Fourier coefficients, r n . Set g(u) = If r(u>) and call its 
Fourier coefficients g n . Then 

£ \n\r n g* n <0, 

n= — °° 

unless r(oj) is constant. 

Proof. Using (49) for the product Rr and the identity immediately fol- 
lowing it, we calculate 



trRT - (2N+ 1) = (2N + 1) [~ f * *y*}r K N fau')duda' - 1 , 

L27T J-Trr(oj') J 



(56) 



where K N (u,w') is the well-known Fejer kernel 3 



sin^ — - — (a) — a) ) 
sir - (o> — <ff) 



It has the following property. Let X(a>) G L 2 (—rr,ir), and let 
X(u>) = C* KNia^XMdu'. 



(58) 



t The kxk matrices A and D have the same eigenvalues. 

X The restrictions in Theorem III and Theorem IV to only a finite number of nonvan- 
ishing components of r(a>) is removed in Section V. 
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Call the Fourier coefficients of X(w) and X(a>), X n and X n , respectively. 
Then 

^n = X„[l-^L_], | n |^2N + l (59) 

= otherwise. 

k 
Thus,ifr(oj) = £ r n e inu , 
-k 

C* K N (<a,u')r(o)du = /■(«') - £ r„|n|e' w (60) 

•/— jr ZiV + 1 — fc 

if only (2N + l)^k. Substituting (60) into the right member of (56), we 
obtain 

k 1 /-»7r pinw' k 

2Zr„|rz| -— I — da'=Zr n \n\g* n . (61) 

-fe Z7T J-wr(w ) -k 

Noting that we have already established that the left-hand side of (56) 
is strictly positive [r(a>) ^ const], the theorem follows. 

V. EXAMPLE AND FURTHER COMMENTS 

A particular example is provided by choosing 

r(a>) = 1 + a cos oj, \a\ < 1, (62) 

and, thus, k = 1. We calculate 



2|_ VT^2] 



2N-l r a -i 

B(N) ((N)= -%==y± + p\>Q, N>0, (63) 



vT 



-1 + Vl-a 2 
P = — 



a 
When N = 0, we have, from (49), 



X max (iV=0) = -- I r(a>)dcoX- f -H*L = — ±— , 
2t J-,r 27rJ-7rr(a)) v 1 — a 2 

and otherwise 

X max (JV)-A + B(iV). (64) 

From (29) it follows that 

sin 2 Q N ^^> (65) 
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where r min and r max denote the minimum and maximum of r(a>), and it 
is interesting to compare numerically (64) with (65). Set a = 0.6, so that 
r mi n/r ma x = 0.25. Then fi N ^ 30 degrees from (65). On the other hand 
AmaxOT = 1.125 + 0.125/9 N , N ^ 0, which means J2n starts at about 63 
degrees and increases to 70 degrees. Equivalently, while (65) allows a 
factor of four between the upper and lower bounds of (4) for this example, 
the more exact evaluation has them differing by only 12 to 25 percent 
depending on NJ 

Exact solutions, as we have just found, may be useful for estimating 
Sl N for some particular r(w). If we already know Qn for some other f (o>) 
and if it is true that there are constants n,n' ^ so that 

-^- ^ Ho>) :g (1 + M')r», (66) 

1 + H 

then, in a similar way to which (65) was derived, we can show that 

^^ < sin 2 fi < (i + M ) (1 + M ')sin2fi N . (67) 

(1 + M )(l + mO 

Equation (67) could be useful when r(a>) has a large or infinite number 
of Fourier coefficients. 

Proceeding further along the direction of bounds, we note that, using 
only r min and r max , (65) can be considerably sharpened. One can show, 
in fact, letting E = r m \Jr mBiX , that 

-***[H(*4)r (68) 

The two basic ingredients are that [see (32) through (34)] 

WA0=maz^ (69) 

and 

JWi^+A-V (70) 

for any positive definite Hermitian matrix. 4 Combining (69) and (70) 
yields 

Xmax(N) ^ SUP |j-TT . (71) 

The right member of (71) is further upper bounded by 

max— f'|M(co)| 2 r(«)rf«X^- f* |w(co)| 2 -—- da>, (72) 
tt(w) 2tt J-tt 2tt J-* r(a») 



t Another bound involving only the ratio r min /r mBX is given in (68). For the present ex- 
ample it yields fi/v ^ 53 degrees, a considerable improvement over (65). 
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where 

— J^*W«)|2d«-l, (73) 

u{o)) being any L 2 function, not just one having Fourier components u n 
restricted to | n | ^ N. 
To maximize (72), consider maximizing 

Q = (Ep.a.) (LP, ^) (74) 

with £p, = 1. Pi = 0, a, > and distinct. Introducing a Lagrange mul- 
tiplier X and differentiating, we obtain 

aiE^- + -i:Pjaj = X (75) 

a, a/ 

for all nonvanishing p/. Whatever values the optimum p,'s take, we may 
regard the sums in (75) as fixed numbers, independent of the index /. 
The resulting quadratic equation in a/ can be satisfied by at most two 
values of a; and, hence, only twop/ are nonvanishing, and are easily seen 
to correspond to the maximum and minimum a/ if we are to maximize 
Q. Also the two pi have equal values. Thus, (for a maximum Q) 



**max — . l fl max > Qmin) ( T I 

^ *^max ^min/ 



(76) 



and (68) follows. 

Our next theorem says something about the limiting behavior of 
^max(N). and in fact bounds the latter away from unity in the general 
case. 

Theorem V. Ifr(a) ^ const. 

lim X ma x(ADl^ri+- L C'rMduxl- f* vVI > 1. (77) 
at— a, 2|_ ztt J-* 2ir J-irr(u>)\ 

This immediately removes the restriction in Theorem III. Also, the 
left side of (56) is now, in general, bounded away from zero, and, by 
simple limiting procedures, Theorem IV also follows without restricting 
r(w) as was previously required. Of course, "— »" is included in the 
statement "<0." 

Proof. We begin with a modified form of (69) (let \p = I» which 
states 

i+pn -p i 
X ma x(N) = max J . T1| • (78) 

Thus, any particular choice for provides a lower bound. We choose for 
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the components <f> k of 0, \k \ ^N,(f>k = &-N,k- Inserting this choice into 
(78) yields 

XmaxW ^ ~ ~ =_ " E gn+Ng-N-mr m -n 

1 -N-N SO n,m=-N 

1 2N 1 T 2JV I 

= - £ rfftr.-t = — E riftr.-i + So 'o • (79) 

gOs,t=0 2g0Ls,t=-2N J 

N 

Since* lim £ g' s g t r s -t = go, the theorem follows. 

-N 
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t Namely, Y£er*-t becomes, in the limit, thesth component of l/r(oj) X r(u) which is, of 
course, 8 h q. 
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